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Abstract. The general form of an integral of motion that is a polynomial 
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1. Introduction 

The purpose of this article is to provide a framework for systematically studying 
finite-dimensional integrable and superintegrable systems with integrals of motion 
that are polynomials of arbitrary finite order, N, in the momenta. In the process, 
we also establish some basic properties of the integrals of the motion and study some 
differences between the integrals in classical and quantum mechanics. 

We restrict ourselves to a two-dimensional real Euclidean plane and to 
Hamiltonians of the form 


H = p\ +pl + V(x,y). 


(1) 


In classical mechanics p\ and P2 are the components of the linear momentum, to which 
we add for use below the angular momentum 


Lz = xp 2 - ypi- 


In quantum mechanics, p\ and P2 (as well as H and L3) will be Hermitian operators 
with 


( 2 ) 


pi = - ihd x , p 2 = -ihdy 


In classical mechanics an Nth order integral of the motion can be written as 


N N-k 




( 3 ) 


k—0 j —0 
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General Nth order integrals of the motion 
or simply 

x = Y1 h’ k ( x ' y)p{P2~ k ~ J > ( 4 ) 

j,k 

with fj t k = 0 for j < 0, k < 0 and k + j > N. The leading terms (of order N) are 
obtained by restricting the summation to k = 0. In quantum mechanics we also take 
the integral of the form © (or ©) but the pi are operators as in © and we must 
symmetrize in order for X to be Hermitian (see Section [3] below). 

We recall that a Hamiltonian with n degrees of freedom in classical mechanics 
is integrable (Liouville integrable) if it allows n integrals of motion (including the 
Hamiltonian) that are well defined functions on phase space, are in involution (Poisson 
commute) and are functionally independent. The system is superintegrable if it 
allows more than n integrals that are functionally independent and commute with 
the Hamiltonian. The system is maximally superintegrable |49j if it allows 2n — 1 
functionally independent, well defined integrals, though at most n of them can be in 
involution. 

In quantum mechanics the definitions are similar. The integrals are operators in 
the enveloping algebra of the Heisenberg algebra H n ~ {aq,..., x n ,pi,... ,p n , h} that 
are either polynomials or convergent power series. A set of integrals is algebraically 
independent if no Jordan polynomial (formed using only anti-commutators) in the 
operators vanishes. 

The best known superintegrable systems are the Kepler-Coulomb system mm 
with the potential V = a/r and the harmonic oscillator ©asm with V = ar 2 . 
According to Bertrand’s theorem [5], the only rotationally invariant potentials in 
which all bounded trajectories are closed are precisely these two potentials. A theorem 
proven by Nekhoroshev m states that, away from singular points, if a Hamiltonian 
system in n dimensions is maximally superintegrable (2 n — 1 integrals of motion) then 
all bounded trajectories are periodic. It follows that there are no other rotationally 
invariant maximally superintegrable systems in a Euclidean space E n . 

A systematic study of other superintegrable systems started with the construction 
of all quadratically superintegrable systems in E 2 and E$ [HEM© Superintegrable 
systems with second-order integrals of motion are by now well understood both in 
spaces of constant curvature and in more general spaces [SHH [2ZH3H El EH EH] • 
Second-order superintegrability is related to multi-separability in the Hamilton-Jacobi 
equations and the Schrodinger equation. The superintegrable potentials are the same 
in classical and quantum mechanics. When commuted amongst each other, the 
integrals of motion form quadratic algebras . 

Superintegrable systems involving one third-order and one lower order integral of 
motion in E 2 have been studied more recently The connection with 

multiple separation of variables is lost. The quantum potentials are not necessarily 
the same as the classical ones, i.e. quantum potentials that depend on the Planck 
constant appear. Some of the quantum potentials obtained involve elliptic functions 
or Painleve transcendents. The integrals of motion form polynomial algebras and these 
can be used to calculate energy spectra and wave functions mmm- A relation with 
supersymmetry has been established [42ll44U46ll56] . It was recently shown that infinite 
families of two-dimensional superintegrable systems exist with integrals of arbitrary 
order 26,33,36,44,53,54,60,61 . Superintegrable systems not allowing separation of 
variables have been constructed E3EE1ES] 
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The present article is to be viewed in the context of a systematic study of 
integrable and superintegrable systems with integrals that are polynomials in the 
momenta, especially for those of degree higher than two. Here we concentrate on the 
properties of one integral of order N in two-dimensional Euclidean space. 

The remainder of the article is organized as follows. Section 2 is devoted to N- 
th order integrals of motion in classical mechanics and includes a derivation of the 
classical determining equations. The determining equations for quantum integrals are 
derived in Section 3. They are shown to have the same form as the classical ones, 
up to quantum corrections. These corrections are shown to be polynomials in the 
square of the Planck constant h and are presented explicitly. The general formulas are 
specialized to low order cases N= 2, ... ,5 in Section 4. Different possible quantization 
procedures are compared and special cases are considered in Section 5. Conclusions 
and a future outlook are presented in Section 6. 

2. Classical Nth order integrals of the motion 

Let us consider the classical Hamiltonian © and the Nth order integral (Q} where 
fj,k{x, y) are real functions. Since X of (0]) is assumed to be an integral of the motion, 
it must Poisson commute with the Hamiltonian 

{H,X} pb = 0. (5) 

The commutation relation Jo]) leads directly to a simple but powerful theorem. 
Theorem 1 A classical Nth order integral for the Hamiltonian m has the form 
Hr! N-2t 

x = J2 ( 6 ) 

t=o j—o 

where fj t 2 t o,re real functions that are identically 0 for j,£ < 0 or j > N — 2 1, with the 
following properties: 

(i) The functions fj, 2 t and the potential V(x,y) satisfy the determining equations 

0 = 2(d x f j - 1 ,2t+d v f j ,2t) - (ti+l)f j+ i,2t-2d x V+(N-2£+2-j)f ji2t - 2 dyV) .(7) 

(ii) As indicated in all terms in the polynomial X have the same parity. 

(Hi) The leading terms in J6]) (of order N obtained for £ = 0) are polynomials of 
order N in the enveloping algebra of the Euclidean Lie algebra e(2) with basis 
{Pi,P2,L 3 }. 

Proof We calculate the Poisson commutator J5j) for H as in JTJ and X as in J3J 

fir vl n®fj,k j+l N-k-j n dfj,k j N-k-j+1 

{H,X} pb = ^-2— p{ p 2 *-2— p{p 2 

j,k 

+ J -Q^f j , k Pi P 2 + [N - k - j)—fj,kP{P2 ■ ( 8 ) 

The first two terms are of order N — k + 1, the second of order N — k — 1. We shift 
j —> j — 1 in the first term and j —> j + 1 in the third to obtain 

n o ( l,fc , dfo,k\ 0 N-k-j+l 

3ik 

+ (U + + ( N -k- AP 2 ~ k ~ 3 ~ 1 . (9) 
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The terms of order N + 1 in the momenta, in m, correspond to k = 0 and are 


dfj-i,o , dfjfi = 0 
dx dy 


( 10 ) 


Eq. (fTT)l) is the condition for the highest order terms of X to Poisson commute with 
the free Hamiltonian Hq = pf +p|, thus proving the third statement in Theorem 1. 
From the form of ©, we see that even and odd terms in X are independent. This 
proves statement 2 of the theorem. Finally, we shift k —> k + 2 in the second term of 
(0. The coefficient of p{p% _fc_J+1 (after the shift) must vanish independently for all 
(j. k) and we obtained the determining equations © ■ 


Corollary 1 The classical integral 0) can be rewritten as 

L#J N-2i 

V _ \ ^ a t N—m—n rn n , \ ^ \ ^ r j N—j—2£ /i i \ 

X ~ / ^ An— m— n,m,n-^3 V 1 P 2 “I" / v ^ Jj i 2iPiP2 5 (H) 

0<m-\-n<N £=1 j =0 

where AM-m-n,m,n are constants. 


Proof: As noted above, the determining equations for fj i0 given by (flUl) do not depend 
on the potential. The solutions of m are 
N-j j 

ho = EE ( N ;? m m )A N - n - m x N -i- n (- y y- m , (12) 

n —0 m =0 

which give the form of the integral (fTTll . Thus for all N the solutions of © for i — 0 
are known in terms of the (N + 1 )(N + 2)/2 constants figuring in (fTTll . 

Let us add some comments. 


(i) For physical reasons (time reversal invariance) we have assumed that the functions 
fj,k{x,y) € R from the beginning. This is actually no restriction. If X were 
complex, its real and imaginary parts would Poisson commute with H separately 
(for V(x, y) G R) and hence each complex integral would provide two real ones. 

(ii) The number of determining equations © is equal to 


["±i] 

E 


e=o 


{N - 21 + 2) = 


j(N + 3) 2 Nodd 

|(AT + 2)(IV + 4) Neven. 


(13) 


For a given potential V(x, y) the equations are linear first-order partial differential 
equations for the unknowns fj, 2 i{x,y). The number of unknowns is 


E (N-2£ + l) = 

£=0 


m±lM±3} Nodd 

i(AC + 2) 2 Neven. 


(14) 


As is clear from Corollary 1, the determining equations for fj t o can be solved 
without knowledge of the potential and the solutions depend on (TV + 1)(N + 2)/2 
constants. Thus, TV + 1 of the functions namely fj t o, are known in terms of 
(TV + 1)(TV -)- 2)/2 constants. The remaining system is overdetermined and subject 
to further compatibility conditions. If the potential V(x,y ) is not a priori known, 
then the system © becomes nonlinear and V ( x , y) must be determined from the 
compatibility conditions. We present the first set of compatibility conditions as a 
corollary. 
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Corollary 2 If the Hamiltonian m admits X as an integral then the potential 
function V(x,y) satisfies the following linear partial differential equation (PDE) 

N-l 

o = £ ) j l(j + i)fj+i,od x V + (N- j) fj,od y V ]. (15) 

i =0 

Proof: This linear PDE is determined by the compatibility conditions for the I = 1 
set of determining equations, namely 

2 (d x fj- 1,2 + d y f jt2 ) - [(j + 1 )fj+i,od x V + (N - j)fjfid y V] = 0. (16) 

Therefore when X is an integral, the functions fj ,2 exist and satisfy msD and so the 
potential satisfies m This PDE depends only on the constants Aj t k,e in the highest 
order terms of X, (ED- 

(iii) For N odd, the lowest-order determining equations are 

fi,N-iV x + f 0 ,N-iV y = 0. (17) 

In particular, for the N = 3 case, the compatibility conditions of m with the 
determining equations for fj^ lead to nonlinear equations for the potential mm- 


3. Quantum Nth order integrals of the motion 


In this section, we shall present a theorem analogous to the classical one of the previous 
section, but applied instead to quantum systems. That is, we would like to show that 
given an Nth order differential operator X that is formally self-adjoint with respect to 
the Euclidean metric and that commutes with a given Hamiltonian, then the number 
of independent functions and determining equations are equal to that of the classical 
case and the determining equations are the same as in the classical case, up to quantum 
corrections which are polynomial in ti 2 . Specifically, we have the following theorem. 

Theorem 2 A quantum Nth order integral for the Hamiltonian m has the form 
1 L#Jat-2£ 

x = 2 £ {fjW’P(pZ~ n ~ 3 }’ ( 18 ) 

£=0 j=0 

where fj^t are real functions that are identically 0 for j, I < 0 or j > N — 2 1, with the 
following properties: 

(i) The functions fj, 2 i and the potential V(x,y) satisfy the determining equations 

0 = M m , (19) 

with 


21 = 2 (d x fj-Ifii + dyfj,2i) 

— ((j + l)fj+l,2t—2d X V + (N — 2t + 2 — j)fj, 1l-ldyV + h 2 Qj,2l) ) 
where Qjjii is a quantum correction term given by 

Qj, 2t = (ffd x f>j-i,2i + + d\(f>jju-i + dy<f>j t 2e-i) (20) 

1-2 2n+3 

E v /_ pA\ n (3 + m \ (N — 21+ 2n + 4 — j — m\ / om o2n+3—mrr\ r 

h > { m H 2n + 3 -m )( d x d y P)/i+m,2^-2n-4 


n —0 m =0 
2^—1 


E v /_fc2\ L(n-i)/2j fj + m\ (N — 2£ + n + 1 + j — m\ / am p>n-m'tr\± 

l^( n ) ( )( n-m )( d x d y V ) ( Pj+m, 2i-n-l, 


n—1 m =0 


General Nth order integrals of the motion 


6 


where the (j>j t k are defined for k > 0 as 

“)( 


1 2b ~ e ( fc2\b—1 


rN - 


21 + 2 b-j 
2 b — e — a 


y a A 


2b—e—a 

y 


fj+a,2t-2b, ( 21 ) 


6=1 a =0 


with e = 0,1. In particular (fj t o = 0, hence Qj.o = 0 so the I = 0 determining are 
the same as in the classical case. 


(ii) As indicated in m, the symmetrized integral will have terms which are 
differential operators of the same parity. 

(in) The leading terms in 1 1 #1) (of order N obtained for 1 = 0) are polynomials of 
order N in the enveloping algebra of the Euclidean Lie algebra e(2) with basis 

{Pi,P2,L 3 }- 


Notice that the parity constraint on the integral m reduces the number of 
possible functional coefficients by about half. Indeed, if the integral is expanded out 
with the derivatives moved to the left then the integral would be of form 

^ = E ~ h2< t>iat)p{P2~ 2e ~ J ~ *^E ^j^-iPipA 2 ^ 1-0 ■ (22) 

1,3 

Thus, for a general, self-adjoint Nth -order integral that commutes with H, 
approximately half of the coefficient functions depend only on derivatives of the 
functions these are the <f>+it- 1 - In general, the functions <pj,k (1211) are polynomial 
in H 2 . 

We begin the proofs by showing that, modulo lower order integrals of motion, X 
can be taken to be self-adjoint. 


Lemma 1 Given X an Nth order differential operator that commutes with a self- 
adjoint Hamiltonian H, then X can be assumed to be self-adjoint. 

Proof. We can always write X = X+ A + A Ty yl . Using this, we obtain 




and its Hermitian conjugate 

0 fl, 


2 

r X-X* 


,H\ 


which together show that the self-adjoint and skew-adjoint part of the operator X 

V XV Y’t 

must simultaneously commute with H. Since, — commutes with H so will i=—p — 
and hence it is possible to assume, without loss of generality, that the integral is 
already in self-adjoint form. Q.E.D. 

Next, we show that a commutator of the form [/, d 3 x dy +k ~I] can be written as 
a sum of lower order anti-commutators which differ in degree by 2. We assume for 
convenience that all of the functions are smooth but of course this could be reduced 
to require only the number of necessary derivatives, in this case k + 1. 


Lemma 2 Given a real f € C°°(R 2 ), there exist real functions, f m , n G C'°°(R 2 ), that 
satisfy the equality 


k min{y n+l} 

[f,did 1 y +k ~ 3 ]=J2 E {fm,n,dt m d k - n - j+m } (23) 

n=0 m=max{0,n—fc+j} 


and furthermore the functions with odd second index, f m , 2 (.+i, are identically 0. 
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Proof. We define the f m ,n inductively on n. Notice that the summation indices on the 
right-hand side of the identity (l2dl) lie within the following regime 

A = {(to, n)|0 < j — to < j, 0 < k — n — j + m < k — j + 1}. (24) 


Furthermore, the index n determines the total degree of the monomial 
d 3 x ~ m dy~ n ~^~ m \ We define the identity operator to be d x dy = Id. 

Consider the quantity 

u = ■'] - £ (25) 

n,mG A 


Expanding out U gives 

U = -2 J2 r '‘.",d'r 
n,mG A 

with 


V ’ 


Un,m — yjm,n ^ (j + n _ m ) °x °y J 

n-l \ 

, _ (j - m + a\ (k - n - j + m + 1 + b - a\ na ol+b-a r ) 

+ a 1 + 6- a )° X °V Jm-a,n-l-b , 

6=0 aGcr ' 


(26) 


where a = {max{0, b — n + m} ... min{TO, k — j — n + b + 1+ to}}. Hence for a hxed 
to, n we have a recurrence relation on the f m ,n s which will make each U mtn = 0 this 
relation is a recurrence on n given by 


Jm,n 


1 

2 



k - j + 1 \ 
1 + n — m) 


d™dl +n - m 


f 


(27) 


1 

2 


EEC 


m + a\(k — n — j + m+ l + b — 
a / \ 1 + 6 — a 


a )d a x dl +b - a f , 


By induction on n, the functions f m ,n exist and are real. The first two functions 
(n = 0) are 

fo,o = -^d y f, fi,o = ~dxf, j> 0, 1 + k — j > 0. 

If j = 0 or 1 + k — j = 0, then either /qo = 0 or /o,o = 0, respectively. The summation 
over a and b ensures that (to — a, n — 1 — 5) (=A and so the functions are well-defined 
recursively on n for all to, n £ A. 

Finally, we see that if k is even then k + 1 will be odd and so the operator 
[/, ftjc+k - -!] w in be self-adjoint which gives the requirements that the sum 


k j 

EE{/^- s i m E _i+ro } 

n=0 m =0 

contains only even terms in k — n and hence only terms with n = 21 can be non-zero. 
Similarly, for k odd, k + 1 is even and so both sides must be skew-adjoint and hence 
k — n must be odd which again gives the requirement that n = 2 1. Q.E.D. 

Next, we use the previous lemma to show that any self-adjoint operator X can 
be put into symmetric form so that the functional coefficients are real. 
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Lemma 3 Given a general, self-adjoint Nth order differential operator, X, there exist 
real functions fj t k such that 

x = ( 2g ) 

k j 

Proof. Given a self-adjoint Nth-order differential operator X, we can always move all 
of the functional coefficients to the left to obtain 

x = (29) 

k j 


Fj'k = 0,Vj<0,k<0,j>k,k>N, 


where the functions Fj^ are possibly complex; we write F J: k = F^ k + iF^ k . Since X 
is self-adjoint, it can be expressed as 

X = i(X + Xt) = i E E d F *k’ d i d v ~ k ~ j } + i\ F lk, did?-*-*]) ( -ih) N ~ k . (30) 

k j 

Next, we use the previous lemma to show that for a given j , k there exists real functions, 
call them gff n such that 


i[F! k dZd?- k -*]{-ih) 


N-k 


E 


j 


{(-l) n h 2n+1 gt k , 2n ,tt~ m P2 


j—m '-N-k—l—2n—j-\-m 


}■ 


n =0 m —0 


For simplicity, we defined gi^ 2 n — 0 whenever (m, 2 n) £ A (PM1) . Thus, there exist real 
functions f h k such that X can be written as 

x = ^EE^h fc >PipE fc_J }- (31) 

k j 


where 


Q.E.D. 


fj,k = 


R 


[(fc-l)/2] j 

E E(-i 


n =0 m—0 


Proof of Theorem 2 First, we prove that, modulo lower-order integrals of the 
motion, X contains only even or odd terms. By Lemma [TJ modulo lower order terms, 
X can be taken to be self-adjoint. Then, by Lemma[3]we know that X can be written 
in the symmetrized form 


X =x e + x° 

lflN-2£ 

xe = o E E (-^ N ~ 2£ {h2,,d^- 2 ^) 


t—0 j =0 
l^] N -2£-l 


x °=oJ2 E 


(32) 


i=o j =o 


where the X e and X° have the opposite parity under time reversal (complex 
conjugation). Thus, since Ft is completely real, we know that X e and X° must 
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independently commute with H and so, modulo the lower-order integral X°, X can 
be written as 

1 l%]N-2e 

x = ~ E E (-in) N ~ 2e {h2e, did?- 2 *-*}. (33) 

e=o j =o 

Next, we consider the determining equations for the functions fj,k and the 
potential V making X a quantum integral. Take A' and H as above both self-adjoint 
and write 

l+N 1+N-k 

[X,H}=J2 E M j>k didl +N - k -^-ihf +N - k . (34) 

k —0 j =0 

We will show that the terms Mj, 2 t+ i can be written as differential consequences of the 
Mj t 2 e and that the requirements 0 = Mj, 2 i are given by (flTHl . which are the determining 
equations for the system. We begin by showing that the terms Mj^e+i can be written 
as differential consequences of the Mj, 2 i- We know that [A, H] is skew-adjoint and so 
we have the requirement that 
l+IV 1+iV-fe 

E E (■ M jtk didl +N - k - j {-ih) 1+N - k y + Mj, k didl +N - k - j (-ih) 1+N ~ k = 0. (35) 

k —0 j =0 

Now, from (13311 we know that (— i) N X will be even with respect to complex conjugation 
(time reversal) and so we can see that [(— i N )X, H] will be real. Thus, since we have 
the equality 

1 +N 1+N-k 

l(-i N )X,H] = Y, E () 2+N ~\ 

k —0 j=0 

we see that both side will be real differential operators which implies that the even 
terms, Mj, 2 t, are completely imaginary and the odd terms, Mj, 2 Z+ 1 , are real. Hence, 
we can compute © 

l+N 1+N-k 

0 = E E (Mj,kdidl +N - k - j {-ih) 2+N ~ k y + M jtk dldl +N - k - j {-ihf +N - k 

k —0 j =0 

= E[M jlM , + {Mj^e+i, (36) 

3,1 


The coefficient of d^dy 2£ J '(— ih) 1+N 2£ in (l36l) gives 

0 = Mj t 2t+i 

a. 2 Tt ~\~i / t2\Ti 

, \ n ) (j + m\(N-2(. + 2n + l-j + m\xmp,2n+l-mn,f 

+ ImJ { 2n + l-m d V M j+m , 2 e-2n 


n —0 m =0 

i-1 2n+2 


2 V ■ 
(—h 2 ) n 


fc 2 1 n > fj + m\fN-2e + 2n + 2-j + m\ a ma 2 n+ 2 -mnf 

{ 2n + 2 — m J d * d V M j+m ,2i-2 n - 1- 


n—0 m =0 


For the case, t = 0 we obtain 

—Z /z 

M jt i = — ((j + l)d x M j+ i, 0 + (N-j + 1 )d y Mj, 0 ) 

and so by induction on £, we can see that the odd terms Mj, 2 t +1 are linear combinations 
of derivatives of the Mj, 2 t. 
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The quantity can then be directly computed using the expansion of X, as in 
(12211 with the functions given as in (EH), to obtain HU). Notice that this equation 
agrees with 0 modulo terms which are polynomial in h 2 and which vanish in the 
classical limit (h 0). 

Finally, as discussed above, if the requirements Mj^e = 0 and their compatibility 
equations are satisfied then so too will be Mj^e+i = 0 and so the requirement 
[H. X] = 0 will also be satisfied. 

To finish the proof of Theorem 2, we show that the highest order terms of A' are 
in the enveloping algebra of the Lie algebra generated by p\ , f >2 and L 3 = ypi — xp 2 • 
The determining equations for the f J: o are the same as in the classical case m and 
hence their solutions are the same m- Thus, if such a symmetry operator exists, then 
the highest order terms agree with those of an Nth-order operator in the enveloping 
algebra of the Lie algebra generated by pi,f >2 and L 3 and so it is always possible to 
express the highest order terms as operators in the enveloping algebra, with suitable 
modification of the lower order terms. Q. E. D. 


Just as in the classical case, the highest-order determining equations m can be 
solved directly and the functions fj,o are given by 03- However, as will be discussed 
in Section 4 below, the choice of symmetrization of the leading order terms will lead 
to ^-dependent correction terms in the lower-order functions. The quantum version 
of Corollary 2 still holds and, remarkably, the linear compatibility condition is the 
same. This does not imply that the quantum and classical potentials are necessarily 
the same since further (nonlinear) compatibility conditions exist. They are in general 
different in the two cases. 


Corollary 3 If the quantum Hamiltonian m admits X as an integral then the 
potential function V ( x , y ) satisfies the same linear PDE as in the classical case, namely 

(W- 

Proof: The I = 1 set of determining equations are 

M jt 2 = 0. (37) 

with 


Mj ,2 — 2 (d x fj- 1,2 + dyfj, 2 ) — [0 + l)fj+i,odxV + ( N — j)fj,od y V + h 2 Qj.2] , 
with quantum correction term 

Qj, 2 = ^dx 4 >j-i ,2 + 2 d y (f>j^2 + d 2 (f>j t i + 

The functions , coming from expanding out the highest order terms are 


PjA 


3 + 1 


dxfj+1,0 


N-j 


dyfj, 0 ) 


Pj, 2 


2 

E 

a=0 


j + a 
a 


N - j -a 
2 — a 


d x d l a fj+a, o- 


The linear compatibility condition of (|37[) is obtained from 
N-l 

°= E 

j =0 

N-l 

= E ) j [{3 + l)d x Vf j+1 , 0 + (N- j)d y Vfj, 0 + h 2 Qj' 2 ] ■ 

3 =0 


( 38 ) 
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The coefficient of H 2 in this equation vanish, as the relevant terms are zero, i.e. 

d*- l - J dlQ h 2 = 0, 

since the functions fj .o are polynomials of total degree at most N, see (fT^l) . Thus, the 
potential satisfies the linear compatibility condition (H3, as in the classical case. 

In fact, this result can be obtain directly by considering instead the following 
form of the integral. Consider a general, homogeneous polynomial of degree N in the 
operators pi , P2 and £3 that is self-adjoint, call it Pn{pi,P2, L 3 ). By Lemma 2 above, 
this operator can be expanded out as 

1 . [JV/2J 1 

Pn(Pi,P 2,L 3 ) =-{fj,o,pip%~ J }+ 

t—1 

Thus, the operator X given in m can instead be expressed as 

UV/2J 

X = P N (p 1 ,p 2 ,L 3 )+ ^ ( 39 ) 

t=\ 


In practice Pn(pi,P 2 , M) is often chosen as 


P N (P1,P2,L 3 ) = - ^A N -m-n,m,n{PTP2,^- m - n }, 


(40) 


although this choice is not necessary for what follows. Different possible choices of 
symmetrization will be discussed in Section 4. Let us now consider the determining 
equations for the fj.'ie- Recall the (j. 21) determining equation is obtained from the 
coefficient of p\p 2 ~ 2i+l ~ : ’ in [X, H], Expanding [X,H] gives 


[X,H] = [P n ,H 0 )+[P n ,V]+ 


LJV/2J N— 2 t 

E z. 

t =1 3 =0 




.(41) 


By definition, the first commutator is 0. The third commutator will give exactly the 
determining equations (flTfl) with replaced by fj^t except that there will be no 
fj : 0 terms in the <j)j,k’s. To be precise, expanding out all the terms of X except Pn 
gives 


[A 72 J N -21 


[_N/2\ N-21. 


E E 2 2f J }= E E (/j:2^ - p'p 


N-21 -3 


e=i j=o 


where 


i=i j =0 

LJV/2J N-21-1 


E E hw-ip{P2 2t+1 \ 


(42) 


•6=1 3=0 


21-1 




a =0 


Therefore, the coefficient of fy[p 2 ~ 2t ~^ +1 coming from the third commutator in (1411) 
are —ihMj^t lUTJh with _/£*, replaced by fj.k and 4>j t k replaced by call these 
Mj 2 i- To finish the determining equations, we would need to compute the coefficient 
°f PiP 2 2e+1 ~' i in [-P/v, V], These equations would clearly depend on the choice 
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of symmetrization, discussed later, and even in the case m are not particularly 
illuminating. 

More interesting to consider is the dependence on h of the determining equations 
for Notice that the highest order terms in h are absent from the <pj. 2 t.-Y and 

so, whereas in general the functions <pk, 2 t-e have terms of order £ in h 2 , the functions 
<pk, 2 t-e will have terms of only £ — 1. Therefore, while the functions Mj^i and hence 
the determining equations for fj^e have at most order £ in h 2 , the functions Mj, 2 i 
depend at most £—1. As mentioned above the determining equations for fj. 2 t coming 
from the coefficient ofpjp^ J will be —ihMj. 2 t plus the coefficients coming from 

[Pn , V]. These coefficients will be at most order £ — 1 and so will not add to the 
degree of h 2 in the determining equations for the adjusted functions which will 
be at most order £ — 1. In particular the second set of equations, those for fj. 2 , will 
have no quantum correction terms at all and will be the same as in the classical case. 
This is in agreement with the previous corollary since the compatibility conditions are 
independent of the functions fj .2 or fj. 2 - 

4. Special Cases: N from 2 to 5 

4.1. Case N = 2 

This is the most well-known and researched case because of its connection with 
separation of variables. However, it proves illustrative for future examples. After 
expanding, the integral of motion (1151) becomes 
2 1 

X = fj, 0 PlP 2 ~ J - + fo ,2 ~~ ^ 00 , 2 , ( 43 ) 


00,1 — dyfo,0 + 7J<9z/l,0, 01,1 — dxf2,0 + yf 1,0 

1 2 

0°,2 = ^^0“0 2 a fa, 0, 

a =0 

in agreement with Ell) . Notice that these all depend on the higliest-order terms o 
which lie in the enveloping algebra of E 2 and are given by m- From these solutions 
we can immediately see that any third-order derivative of a function / J: o is identically 
zero. 

The final determining equations are Mj t 2 = 0 with 

Mj, 2 = 2 dxfj-^2 + 2 dyf jt 2 - ( j + l)f j + i,od x V - (2 - j)fj, 0 d y V 

- K 2 (2d x <f)j-i t 2 + ^dy<t>j,2 + d 2 0j, 1 + dy(j) jt 1 ). (44) 

Note that the quantum correction term (the second line) depends on third derivatives 
of the fjfl and hence vanish on solutions m ■ Thus, the determining equations for 
second-order integrals of the motion are equivalent in both the classical and quantum 
cases. They reduce to 

Afo,2 = 2d y fo,2 — fi.odxY — 2fo,od y V = 0 
Afi ,2 = 2d x fo t 2 — ‘Zf2fid x V — fifidyV = 0. 


(45) 

(46) 
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4.2. Case N = 3 


We now turn to the case N = 3. This case was investigated in [T71IT8] . where the 
determining equations were given. Here we obtain the same results, presented in the 
notation of this paper. 

As for all N, the (j>j.o are identically 0. There are thus essentially two families of 
the (f>'s : those that depend only on /y 0 


,d = £ 


a =0 
2 


* = £ 


a —0 


j + a 
a 

j + a 
a 


3 - j ~a 

1 — a 

3 - j ~a 

2 — a 


d x d v a fj+a ,0 


d x d y a fj+a, 0 - 


and one that also depends on fj t 2 , 


2 2b~ 1 /_fc 2 \ 6—1 

<(>0,3 = £ £ 4-26- 


6—1 a —0 


(47) 

(48) 


(49) 


Let us now consider the quantum correction terms. As always we have Qj t o = 0 (1211) . 
The determining equations M, j3 = 0 are given by m with solutions m- The next 
set of quantum corrections are 


Qj, 2 — (29 x ^j_i ! 2 + yi + . (50) 

While it is not immediately clear that these vanish on solutions of m, one need only 
observe that each of these terms involve third-order derivatives of the fj i 0 all of which 
vanish except for 


^z/o.o — — dxdy f 1,0 — d x dyf2,o — —dyf3,0 — 6 A 30 o. 
From these, it is straightforward to compute the only non-trivial terms 


d x <t>j~ 1,2 — dx&j,! — 


0 j = 0 

— 6 A 300 j = 1 


(51) 

(52) 


dy&j ,2 ~ dyftjA 


0 j = 0 

6 A 300 j = 1 


(53) 


Thus, Qj t 2 = 0 on solutions of m and so the next set of determining equations are 
again independent of h and given by 

d y fo ,2 = ^fi,od x V + ^fofidyV (54) 

d x fo,2 + d y fi,2 = f2,od x V + fifidyV (55) 

d x fi,2 =\f3,AV+ l -f 2 ,od v V. (56) 

To compare these equations to those of Gravel El, the potential is multiplied by a 
factor of 2 (to make up for the choice of a factor of 1 in the kinetic energy term in 
HI)) and the functions f jt0 are related to f k via / 3 , 0 = / 1 , / 2 ,o = / 2 , /i,o = /3 and 
/o,o = f±- The other set of functions are related via g\ — /i ,2 — 5H 2 yA 300 + |fi 2 A 2 io 
and 52 = fo .2 + 5fi 2 a;A3oo + |S 2 A 2 oi- The quantum correction terms are related to the 
choice of symmetrization, as discussed in Section 5. 
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Turning now to the final determining equation. The quantum correction term is 

Qoa = [p x 4> o>3 + 9y(f> 0 , 3 ) 

321 

- E (d*dy~ m V)fm ,0 ^ E Pxd 2 v - m V)ct> mtl - E 

m=0 m —0 m =0 

Note that from this expression it appears that Qo ,4 has a term depending on h? (from 
(^ 0 , 3 ) which would lead to an term to the final determining equation Mq 4 . However, 
if we inspect this term, we can see that it contains only fifth-order derivatives of the 
functions /yo and so will vanish, see m- This leads to the simplification 

0 = - fi,2d x V — fo,2dyV 

- & ( - \ E (9Td 3 y- m V)f m , 0 - l&dyhMV - \ (d x dyf lfi ) dyV. j (57) 
V m=0 J 

Note that this equation is linear and homogeneous in derivative of the potential V. 

Summarizing the results from this section, there are three families of determining 
equations. The first ensures that the leading order terms are in the enveloping algebra. 
This set as well as the second set of equations are the same in the classical and quantum 
case. The final set of equations (in this case one equation) does have a quantum 
correction term, linear in h 2 which is also linear and homogeneous in the derivatives 
of the potential V. 


4.3. Case N = 4 


The structure of the fourth- order integrals is similar to that for N = 3, except that the 
number of equations in each set will be different. These determining equations were 
obtained in Ref. 1551. for completeness we present them here. There are 4 equations 


for the fj t 2 functions. 

= -^fi, od x V + 2f 0 , 0 d y V + H 2 (6yA 400 — -A 31 o) 


dyfo,2 

(58) 

dxfo,2 

+ &yf\,2 = f2fid x V + —flfidyV -f- (6xA 4 qq + —A 3 oi) 

(59) 

dxfl,2 

+ d y f 2 ,2 = -^fsfidxV + f 2 ,od y V + h 2 (6yA 400 — -A 31 o) 

(60) 

d X f 2 ,2 

= 2f4,od x V + —f 3 fid y V + H 2 (6 xA 4 qo + — A 3 qi). 

(61) 


Notice that, unlike the case N = 3, the quantum corrections Qj , 2 are not identically 
0. However, the equivalent computations obtained in [55] do not have any quantum 
corrections. As described above, this is due to the fact that the integral in m was 
not assumed to be of the form m but instead as in 


A = P 4 {p 4l p 2l L 3 ) + i{/ 2 , 2 ,P?} + ^ {/l ,2 j PlP 2 } + 7 }{fo, 2 ,pl} + /o,4, 
with 


3 1 

fo ,2 = fo ,2 — h 2 ( 3 y 2 A 4 oo — -A 310 y + -cia : 2 C2X ° 3 


(62) 

(63) 


fi ,2 = fi ,2 ~ ft 2 ( 6yxA 400 - axy - 3zA 3 i 0 - 2c 4 a; + ^yA 30 1 - c 2 y + c 5 


(64) 


/ 2 ,2 — / 2,2 — P 2 ( 3 X 2 A. 


3 . 1 2 

*400 + ^A 3 oi + -Cl?/ + C ±y + c 6 


(65) 
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where the coefficients Cj depend on the choice of symmetrization. 

For the final set of determining equations, there are two that need to be satisfied, 
instead of one as in the case N = 3. They are given by 

0 = 2d x f 0A - (2h a d x V + f lt2 d v V + h 2 Q 1A ) , (66) 

0 = 2d y fo A - (fi,2dxV + 2f 0 , 2 dyV + H 2 Q 0A ) , (67) 

with 

Q 1,4 = Fyyy/1,0 214yy/2,0 314xy,/3,0 414xx,/*4, 0 

( 9 3 

— Vyy ( ~(6® 3 + 12x2/)Aioo H— GxyA ^ 10 — (77X 2 — 3 y 2 )Asoi — 2 yA 2 u — 3xA.202 + 2xA.22o — 77^103 + ^121 


— 2V xy ( (12 x 2 y — 6y 3 )A 400 — (3x 2 — -y 2 )A 3 io + QxyA^i + 2yA 2 Q 2 — 2>yA 22 Q — 2 x^ 2 n + -A 130 — ^-112 


- 314 : 


3 1 

—6XJ/A400 + 3xy^43io — —A301 — xA 2 q 2 + yA 2 n — —A\ 2 i 


— V y ( 48 x 2/^400 + 122 /A 301 — 12x^-310 — 4A.2 h) 

— 214 ((6x 2 — 182/ 2 )^44 oo + 92/^310 + 3x^301 + A 2 q 2 — 3 A 220 ) 

+ dyd x fo ,2 + '^d y d 2 Ji , 2 + 2 d 3 J 2 ,2 + 7 ^/ 1,2 + d 2 d x f 2 , 2 

Q 0,4 = AVyyyfoQ 314yy/l,0 2\f XX yf 2 ,Q V XXX f3,0 

3 


1 


— 31y y I —6x 2/^400 + 77 ® A 310 — 3 x 2/^301 + ®A. 2 n — 2/^202 + 77-^112 


— 214 


xy 


(6x 3 + \ 2 xy 2 )AiQo — (-x 2 _ 3 j/ 2 )^43oi — 6X2/A310 — 3x^202 + 2x^220 — 2 yA 2 n — -A403 


1121 


— 14a: ^(12x 2 y — 62 / 3 )^ 4 oo + QxyA^Qi — (3x 2 — t ^ 2 )^ 310 + 2yA 2 Q 2 — 3yA 22 Q — 2 xA. 2 ii + -A 130 — dn 2 
+ 214 ((18x 2 — Qy 2 )A/iQQ + 9x^301 + SyAziQ + 3 A 202 — 4 I 220 ) 

— 14 ( 48 x?/Aioo — 12 x ^310 + I22/A301 — 4^211) 

+ 2d 3 fo, 2 + -d 2 d x fi, 2 + d y d 2 f 2 , 2 + d y d 2 fo, 2 + -<9 3 /i,2 

Again, these determining equations could have quantum corrections quartic in h but 
instead depend only on h 2 . To compare these determining equations with those of [ 551 . 
we use the differential consequences of (15811611) . 


44. Case N = 5 

The integral we would like to investigate is 

2 5 — 2 £ 

x = \ E E ditf- 2 *-*}- (es) 

j =0 

The highest-order determining equations are the same as in the all cases m ■ with 
solutions (USD- However, unlike previous cases, the next family of determining 
equations Mj, 2 = 0 contain non-trivial quantum correction terms, given by 

Qj ,2 = (2d x f>j-i, 2 + 2d y (j)j, 2 + d 2 <f>j, 1 + d 2 (f>j A ) . 


(69) 



General Nth order integrals of the motion 


16 


Unlike the lower-order cases, these quantum correction terms do not vanish and the 
determining equations are 


dyfo,2 

dxfo,2 + dyfl,2 = 


■zh'OdxV + T;fo,od y V - h 2 -A 311 + §A iw x - 3QA b00 xy - QA 4m y 


f2,od x V + 2fi } od y V 

— h 2 (30 A 3 Qoy 2 — 30 Asoo ^ 2 — 12 Aiioy — 12 A 4 qiX + 3 H 320 — 3 A 302 ) 


dxfl;2 + dyf2,2 ~ -;f 3 ,od x V + -f- 2 ,odyV, 
d x f2,2 + dyf3,2 = 2/4,(AU + fzfidyV 

— h 2 (30 A 3 ooy 2 — 30 21 , 500 a ; 2 — 12 A^y — 12 A 401 X + 3 A 3 2 o ~ 3 ^ 302 ) 


dxh;. 


5 1/ 3 

= 7 ;f 5 ,od x V + -U,od y V - h 2 ( 6 A 401 y + 30 A 500 xy - -A 3 u ~ 6 ^4io£ 


Note that while these equations contain possibly non-trivial quantum corrections, they 
can be removed by taking the form of X as in (1391) . We omit the remainder of the 
equations as they quickly become unmanageable, although it is interesting to note 
that in both Mj t 4 and Afo ,6 the highest degree of h is 2 and 4 respectively; that is, as 
in the cases for IV < 3, the highest-order terms in h are missing from the two lowest 
families of equations. 


5. Quantizing Classical Operators 

In this section, we would like to make a few comments and observations about the 
relationship between classical and quantum integrals. It is clear from the previous 
sections that, given a quantum integral of motion of the form in (1181) that commutes 
with a Hamiltonian H , and assuming that both the potential V and the functions 
fj_ 2 i are independent of h, then the classical integral 

[f ] N—2& 

X = fx 2 UV > iP 2 ~‘ ll ~\ ( 7 °) 

(.—0 j =0 

will be an integral of the motion for a classical Hamiltonian with the same potential. 

However, the implication is clearly not reversible. This observation and the 
general question of quantizing a classical system has received much interest over 
the years, see for example [22H24] . One of the first questions that appear in the 
transition from the classical to the quantum case is the choice of symmetrization, 
assuming cannonical quantization. Clearly, the symmetrization used in (1181) is not 
unique. Indeed, there are many possible choices. The most general choice would be 

fj,MPiP2~ 2e ~ J S(fj, 2 t,c a , b ) (71) 

with 

N-21 j 

Of £ \ _ \ ^ \ ^ (-a '■b—a £ a ~N-2t—j—b-\-a . - 7 —a ~N-2£—j—b-\-a r -a ~b— 

S{fj,2e,C a , b ) = 2 ^ 2^ Ca ’ b \ P lP2 fiWPl P2 +Pl P2 Jj,2tPlP2 

b— 0 a—0 

with Y^ab Ca h = If is true, however, that the choice of symmetrization does not 
affect the general form of the integral and furthermore the different choices will lead 
to quantum correction terms in the f h -u as polynomials in h 2 . In particular, given a 
choice of symmetrization as in CD. we shall show that this choice is equivalent to 
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the standard one chosen in (fT%l) with quantum corrections to the coefficient functions. 
The determining equations for any choice of symmetrization in (ED are equivalent 
to the standard one chosen in (TTSl) up to appropriate modifications in the quantum 
corrections Qj^e- The formulas given in Theorem 2 only hold for for the canonical 
choice dH). 

Theorem 3 Let fj,u G C°°(K 2 ) be polynomial in H 2 , then the self-adjoint differential 
operator S(fj^e,c a ,b) defined in J7i| j can be expressed as 

[N-2t\ 

S(fj,2t, C a ,b ) = ^{fj,2Z,riv2~ n ~ 3 } - & 

k—0 

where the g^ 2 i+ 2 k are polynomial in h 2 . 

By Lemma [3l the differential operator S(fj^e,c a ,b) can be expressed as the sum 

Qi,b) = — (73) 

m n 

Because of the requirement c a ,b = the leading order terms will agree, g : . 2 f. = fj, 2 i- 
The next order term will have a different parity under time reversal (complex 
conjugation) and so will vanish identically leaving a differential operator of degree 
2 less multiplied by a factor of —h 2 . Thus, by induction, we see that ED holds. 

Note that the case of quantizing classical integrals would correspond to the 
functions fj, 2 t being independent of H 2 , in the theorem. In what follows we give 
an example of a classical system where the choice of symmetrization is slightly non- 
intuitive. 


r ~N—2l—2k—j't \ /'7Q\ 

{ gj,2t+2k,PlP2 }), (72) 


5.1. Example I 


Let us consider an example of a potential that allows separation of variables in polar 
coordinate [52], has a nonzero ^ 2 (^ 3 ) 2 term, and has a non-zero classical limit. This 
example is in fact second order superintegrable and the third order integral given 
below can be obtained from the lower order integrals. The potential is 

a a i j_ “ 22 / 

r 2 


V = 


% + . ,_ 

x 2 x 2 \/ x 2 + y 2 


(74) 


\/x 2 + y 2 

and it allows separation of variables in polar and parabolic coordinates. The classical 
third-order integral is given by 

X = P2L3 + fl,2Pl + fo,2P2 


(75) 


with 


/i,2 = - 


/o ,2 = 


x(ay - a 2 ) 

2\Jx 2 + y 2 

r 2 + 2(322/ 


OL2V 


2y/a 


2/ 


x 2 \Jx 2 


+ 


ai(x 2 + y 2 ) 


The correct symmetrization which keeps /j 2 and /o ,2 fixed is given by 

X = -(P2L5 + 2L3P2L3 + L5P2) + ^ {/o,2 5 P2 } ■ (76) 
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This integral can be expressed in the standard form m via 

/o,o = z 2 , fi,o = -2xy, f 2 ,o = y 2 , /3,0 = o 

and a quantum correction to /o ,2 of 7/4 H 2 so that (17(1 becomes 

3 _ 

X =Y1 + ^{fl, 2 ,Pl} + ~,{fo ,2 + jh 2 ,P 2 }- (77) 

1=0 

Note that the appropriate symmetrization (keeping the lower order-terms free of h 
dependent terms) is neither IT 2 , vo\ as in |17U62I nor simnlv i ft n. dld'X^ 1 as in m- 


5.2. Example II 


In this section, we consider quantum Hamiltonian systems associated with exceptional 
Hermite polynomials of Marquette and Quesne [32>]. We mention also that these 
systems are composed of ID exactly solvable Hamiltonians for exceptional Hermite 
polynomials discovered and analyzed by Gomez-Ullate, Grandati, and Milson [14] . 
The first example of a rational extension of a simple Harmonic oscillator is given 

by 

tt i-2 / n2 o2\ 2/ 2 2\ Sh Uj( 2 iOX ft) 

H = -h 2 (d 2 + d 2 y ) + u 2 (x 2 +y 2 ) + ' 


(2 cox 2 + ft) 2 ' (?8) 

The Hamiltonian admits separation of variables in Cartesian coordinates and hence a 
second-order integral of the motion 

v i 2 2 

Ai =p 2 + u) y . 

Additionally, there are two, third-order integral of the motion given by 


-A 2 = ^{L 3: pip 2 } + {uj 2 x 2 ,L 3 } 

(, I2hu>y(2ujx 2 — H) „ r 12hujy(2u)x 2 — h) „ , 

+ h {2 “V + n ,Pi} - > 2 } 


(79) 

(80) 


x 3 = Li 


3ft 

2ui 


(2 tux 
:{L 3 ,X 1 } 


w 


+ h ( 2L 3 — { 


8huj 2 y 3 (2ujx 2 + h) 


(2o 


Pi} + { 




8Huj 2 xy 2 (2ujx 2 + 3 h) 


(81) 


,P2} 


(2wx 2 + h ) 2 ’ (2w;r 2 + ft) 2 

This system can thus be considered as a quantum deformation of the harmonic 
oscillator. As demonstrated by the existence of the third-order integrals, this system 
falls in the classification of Gravel HU- The connection can be directly obtained by 
setting w = —H/(2a) 2 to obtain the potential V e in the classification [TT]. 

Another example is based on the fourth Hermite polynomial. The associated 
superintegrable Hamiltonian is given by 

TT t2 / q2\ 2, 2 2\ 768 h A 0J 2 X 2 16ft 2 w (2uix 2 - 3ft) 

vl (4w 2 a: 4 + 12wa; 2 ft + 3ft 2 ) 2 4w¥ + 12wa; 2 ft + 3ft 2 

The Hamiltonian admits Aft, as in GD, as well as the following third-order integral 
*2 = {L 3 ,p\} + {u 2 x 2 , L 3 } 


6luL 3 — 


f 24hioy(8uj 2 x 6 + 12ftura; 4 + 18ft 2 wa; 2 — 9ft 3 ) 

\ 4w 2 ir 4 + 12ftwa; 2 + 3ft 2 ) 2 

f 8hujy(8uj 2 x 6 — 12hoj 2 x 4 — 6 H 2 ujx 2 — 27 ft 3 ) 


Pi 


V 


4w 2 ir 4 + 12ftu>a; 2 + 3ft 2 ) 2 


,P 2 


(82) 
















General Nth order integrals of the motion 


19 


Unlike the previous case, this potential is not immediately recognizable in Gravel’s 
classification. A remarkable fact is that this potential is associated with the fourth 
Painleve equation. Indeed, the ^-dependent part of the potential W(x) = V (x, y) — 
oj 2 y 2 + iooH satisfies 

-fi 2 IU (4) - \2u 2 {xW)' + 3(IU 2 )" - 2 uj 2 x 2 W" + 4cu 4 x 2 = 0, 

which is equivalent to Eq. (17) of the paper [17] by scaling the potential and the 
variable. It can also be shown that the potential of the previous Hamiltonian (1781) is 
also a solution of this non-linear equation. Thus, these two systems whose wave 
functions are given by exceptional Hermite polynomials have potentials that can 
be expressed in terms of rational solutions to the fourth Painleve equation [T9|[2(J1. 
Of course, many such particular solutions to the Painleve equations exist but their 
connection to exceptional orthogonal polynomials as well as the harmonic oscillator is 
quite remarkable and will be investigated in future work. 

6. Conclusions 

The main results of this article are summed up in Theorems 1 and 2. They present 
the determining equations for the coefficients of an Nth order integral of the motion 
X in the Euclidean plane E 2 in classical and quantum mechanics, respectively. Both 
the similarities and differences between the two cases are striking. 

The number of determining equations Mj^e = 0 to solve and the number of 
coefficient functions Fj 2 t(x,y) to determine is the same in the classical and quantum 
cases (see m and (HI). If the potential V(x, y ) is known, the determining equations 
are linear. If the potential is not known a priori then in both cases we have a 
coupled system of nonlinear PDE for the potential V (x, y) and the coefficients fj. 2 ( 
(0 < £ < LyJ 0 < j < N — 2£). In both cases the functions / j t 21 are real, if the 
potential is real and the quantum integral is assumed to be a Hermitian operator. 
In both cases the integral X contains only terms of the same parity as the leading 
terms (obtained for £ = 0). The leading terms in the polynomial he in the enveloping 
algebra of the Euclidean Lie algebra. This is a consequence of the fact that the 
determining equations = 0 for £ = 0 do not depend on the potential and are the 

same in both the classical and quantum case. For N > 3 the quantum determining 
equations with £ > 1 have quantum corrections. For £ = 1 the determining equations 
Mj 2 = 0 for j in the interval 0 < j < N must satisfy a compatibility condition. This 
linear compatibility condition (El) for the potential to allow and Nth order integral is 
the same in the classical and quantum case. However, if we consider the determining 
equations 0 and HU respectively for £ > 2 they will differ in the two cases. Moreover, 
new compatibility conditions on the potential arise for each higher value of £. They 
will be nonlinear equations for V (x, y) and will be considerably more complicated in 
the quantum case than in the classical one (see jT71ll8lH0 1148115^115^] for the case N = 3 
and Section 3 of this article for N > 3). 

As mentioned above, the number of determining equations to solve for any given 
N is actually smaller than given by m since those for ff o have been solved for general 
N (there are N + 2 such equations). Even so, solving the determining equations for 
N > 2 is a formidable task, even for N = 3. A much more manageable task is to use the 
determining equations in the context of superintegrability. In a two-dimensional space 
a Hamiltonian system is superintegrable if it allows two integrals of motion X and Y, 
in addition to the Hamiltonian. They satisfy [A, H] = [Y. H] = 0. Assuming that X 
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and Y are polynomials in the momenta and that the system considered is defined on 
E 2 , both will have the form studied in this article. The integrals H, X and Y are 
assumed to be polynomially independent. The integrals X and Y do not commute, 
[X, Y] 0, and hence generate a non-Abelian polynomial algebra of integrals of the 
motion. 

The case that has recently been the subject of much investigation is that when one 
of the integrals is of order one or two and hence the potential will have a specific form 
that allows separation of variables in the Hamilton-Jacobi and Schrodinger equations. 
The potential V ( x , y) will then be written in terms of two functions of one variable 
each, the variables being either Cartesian, polar, parabolic or elliptic coordinates. 
Once such a potential is inserted into the determining equations, they become much 
more manageable. 

For N = 3 and N = 4 pJHH SHI S3 SHI H2 M\ this procedure, at least for 
Cartesian and polar coordinates leads to new superintegrable systems, not obtained 
when both integrals X and Y are of order N < 2. In classical mechanics these 
potentials are expressed in terms of elementary functions or solutions of algebraic 
equations [IBIHBlIflZ] . In quantum mechanics one also obtains “exotic potentials” that 
do not satisfy any linear PDE, i.e. the linear compatibility condition m is solved 
trivially, see |40j . These exotic potentials are expressed in terms of elliptic functions 
or Painleve transcendents [I7|(I5Hi5|[B2]. 

This has so far been done systematically for N = 3 and iV = 4. The formalism 
presented in this article makes it possible to investigate superintegrable separable 
potentials for all N. Another application of the presented formalism is to make 
assumptions about the form of the potential and then look for possible integrals of 
the motion. Hypotheses about integrability or superintegrability of a given potential 
can then be verified (or refuted) by solving a system of linear PDEs. 

Alternative approaches to the construction of superintegrable systems in two 
or more dimensions exist. In quantum mechanics they typically start from a one¬ 
dimensional Hamiltionian H\ = pf + Vj (x) and a Hermitian operator A'i that satisfies 
some predetermined commutation relation with the Hamiltonian. A very recent 
article m was devoted to the “Heisenberg symmetries,” when the commutation 
relation was chosen to be = 1. Combined with a similar relation in the 

y coordinate [H 2 , X 2 J = 1, this leads to interesting superintegrable systems. For 
other assumptions, e.g. that A'i is a one-dimensional ladder operator, see [451146] and 
references therein. 

There are many avenues open for further research. In addition to the already 
discussed applications to Nth order superintegrability already in progress, we mention 
the extension of the theory to spaces of non-zero curvature m and to higher 
dimensions [551f57] . 
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